Nuclear modification of vector boson production in proton-lead collisions at the LHC 



Zhong-Bo Kangi'B and Jian-Wei Qiu^'^.Q 

^Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545, USA 
^Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA 
'^C.N. Yang Institute for Theoretical Physics, Stony Brook University, Stony Brook, NY 11794, USA 

(Dated: January 1, 2013) 

In anticipating the upcoming proton-lead run at the LHC in the near future, we present predic- 
tions for the nuclear modification factor of transverse momentum spectrum of production and 
transverse momentum broadening of vector boson {J/ip, T, W/Z'^) production in proton-lead colli- 
sions at y's = 5 TeV, respectively. We find that the measurement of nuclear modification factor of 
Z° production provides a clean and unambiguous test of the nuclear anti-shadowing proposed in the 
^SJ , recent EPS09. In addition, the dramatic difference in transverse momentum broadening between the 

T-H . heavy quarkonium and W/Z" production could be a signature prediction of QCD parton multiple 

' scattering, and it provides valuable informations on color neutralization of a produced heavy quark 

pair when it transmutes into a physical quarkonium. 
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■ Extensive studies of ultra- relativistic nucleus- nucleus collisions carried out by the experiments at the Relativistic 
I Heavy Ion Collider (RHIC) at BNL and the Large Hadron Collider (LHC) at CERN have demonstrated that the 

hot and dense matter of quarks and gluons, the quark-gluon plasma (QGP) was formed in the interactions of heavy 
^ nuclei at high energies jl), i2i]- In order to extract precise information on QGP properties, one needs to calibrate hard 
. probes by producing them in better-controlled cold nuclear matter. The upcoming proton-lead (p+Pb) collisions at 

■ the LHC should accomplish precisely that, thereby provide crucial benchmarks for understanding the characteristic 
of QGP produced in heavy ion collisions at the LHC [3]. 

I In anticipating this upcoming p-|-Pb run at the LHC in the near future, we propose in this paper two measurements 

to better derive nuclear parton distribution functions (PDFs) and to study the dynamics of parton multiple scattering 
in the cold nuclear matter. In the first part of our paper, we will study the nuclear modification of transverse 
' momentum spectrum for neutral heavy boson production in p-|-Pb collisions comparing with that in p-l-p collisions, 
, and demonstrate that the modification serves as a direct measurement of nuclear modification of gluon distribution. 

■ Then, we study the transverse momentum broadening of vector boson ( J/V', T, W/Z'^) production in p+Pb collisions 
^sj to explore the dynamical consequences of parton multiple scattering in QCD and the formation of heavy quarkonia 

. from produced heavy quark pairs. 

Production of Z° bosons, particularly the full transverse momentum spectrum in p-l-p collisions, has been a test 
ground for short-distance dynamics of QCD and an ideal observable for perturbative QCD treatment of hard probes. 
This is because of the nature of two observed hard scales, the Z^-boson mass Mz and its transverse momentum 
. Pt, and our ability to measure a very wide range of px- When px ^ Mz, Z'^-boson production is a very clean 
r\ • probe of perturbative QCD at a distance scale as small as an attometer. On the other hand, when the two scales 
, are very different, the pT-distribution of Z^ boson production could be a perfect laboratory to study the richness of 
QCD dynamics covered by these two scales as well as the infiuence of non-perturbative physics. At low transverse 
momentum px <C Mz, the pT-distribution, calculated in the conventional fixed-order perturbation theory, receives 
a large logarithm, \n{M^/p'^). One could get as many as two powers of the logarithm for every power of coupling 
constant when the calculation goes beyond the leading order. Therefore, at a sufficiently small px, convergence of 
the conventional perturbative expansion is impaired, and these large logarithms must be resummed. A resummation 
formalism, known as CoUins-Soper-Stcrman (CSS) formalism, was proposed, developed and improved over the years 
^^d has become one of many successful examples in perturbative QCD toolkit. 
The study of nuclear modification of Z^ production has finally become available at the LHC llOl-flal . Nuclear 
modification of the inclusive production and its implication on nuclear PDFs was studied in Ref. [16|. There are 
also studies on nuclear modification for low mass virtual photon production but at high transverse momentum [Tt']. 
In this paper, we will use the CSS resummation formalism to study nuclear modification to the p^-distribution of 
Z'^ production in p+A collisions at the LHC. In proton-nucleus collisions, there are two major nuclear effects that 
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one needs to consider carefully. One is the power corrections coming from the initial-state multiple parton scattering, 
which could be enhanced by the nuclear size and thus potentially could be a large effect. The other is from nuclear 
PDFs within the factorized leading-twist formalism [18l - l2]| . We will show below that the nuclear size enhanced power 
corrections (or the high-twist effect) are small for Z° production at the LHC energies. But, the nuclear modification 
of Z° production is directly tied to the nuclear modification of PDFs, and thus, provides very good constraints to the 
nuclear PDFs. 

One of the sensitive observables to the parton multiple scattering is the transverse momentum broadening 
[23 |. For colorless vector boson production (W/Z^), the incoming (anti-)quark from the proton can undergo multiple 
scattering with the soft partons inside the nuclear matter before the hard collision to produce the vector boson. 
We call them "initial-state parton multiple scattering", which will lead to an accumulative change to the averaged 
transverse momentum square (p'^). This phenomenon is often referred to as transverse momentum broadening. On 
the other hand, heavy quarkonium production in proton-nucleus collisions could have both initial-state and final-state 
multiple scatterings in nuclear medium since the heavy quarkonium is very unlikely to be produced at the moment 
when the heavy quark pair was produced [1^ . The final-state multiple scattering between the produced heavy 
quark pair and the nuclear medium before the pair transmutes itself into a physical quarkonium is sensitive to the 
color, spin and other characteristics of the pair, as well as the hadronization mechanism of heavy quarkonia. That is, 
the nucleus in p+Pb collisions could serve as a femtometer size detector or filter to gain much needed information on 
the formation of heavy quarkonia. The additional final-state interaction for heavy quarkonium production could lead 
to more broadening in quarkonia's transverse momentum. The amount of additional broadening is determined by 
QCD multiple scattering and the formation mechanism of the quarkonia. The dramatic difference in the transverse 
momentum broadening between heavy quarkonium {J/ip, TC") and the WjZ^ production, as presented later in this 
paper, is a signature prediction of QCD calculation of multiple scattering and the role of color. 

The rest of our paper is organized as follows. In Sec. II, we study the nuclear modification of boson production 
in p-|-Pb collisions at the LHC. In Sec. Ill, we present predictions for the transverse momentum broadening of J/?A, 
T, and WjZ^ productions. We conclude our paper in Sec. IV. 



II. NUCLEAR MODIFICATION FACTOR FOR TRANSVERSE MOMENTUM SPECTRUM OF Z° 

PRODUCTION 

In order to cover the full range of transverse momentum of Z^ bosons, we use the CSS resummation formalism [3| 
to calculate the Z^ production in p+p collisions at the LHC energy, 

daA^B-^Z^^X ^ 1 ( dHe'PT-bw(b,Mz,XA,XB) + YipT,Mz,XA,XB), (1) 

dydp^ [2iiY J 

where W term gives dominant contribution when px Mz, and the Y term is perturbatively calculable and given 
in Ref. Q. The Y term allows the formalism to have a smooth transition from the resummed low pt region to a 
region where pt ^ Mz and the fixed order perturbative QCD calculations work well. In Eq. xa = Mz/\/s 
and xb — AIz/^/s with the rapidity y and collision energy y^s. The W in Eq. ([T]) is given by [3] 



w/L 7,,r ^ J W^''''\b,Mz,XA,XB) b<b„,ax /r,^ 

W[t>,Mz,XA,XB)-<, W^crt^f,^^^^Mz,XA,XB)F^Pib,Mz,XA,XB;braa.) b > braa. ^ ' 



where b„iax ~ l/(few GcV) is a parameter to specify the region in which 14^^°'* is perturbatively valid, and 
is a non-perturbative function that determines the large b behavior of W and is defined below. In Eq. ([5]), the 
W^°'^*^{b, Mz, xa,xb ) includes all powers of large perturbative logarithms resummed from ln(l/6^) to ln(Af|) and has 
the following form ^ 

lyP'^^H^ Mz, XA, xb) = e-^^"'*^-) c/b, xa,xb) , (3) 

where c is a constant of order one [1, 0], and S{b, Mz) — In ("Tlf ) ^(q^s(/^)) + B{asip)) , with perturba- 

tively calculated coefficients A{as) and B{as) given in Ref. 7] and references therein. In Eq. ([3]), M^^°''*(6, c/b, xa,xb) 
has no large logarithms and is given by 

W^"\b,^,XA,XB) = (Tq fi/A{xA,lJ-^J^)h/BixB,fJ-=^), (4) 

i=q,q 
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where (Tq is the lowest order partonic cross section for a quark and an antiquark to produce a boson and the 
functions fi^j^ (and fj^s) ^^'^ the modified parton distributions given by 0] 



„ JXA >• > 



(5) 



where runs over all quark flavors as well as gluon, (pa/Ai^^ m) ^'^^ the normal proton or effective nuclear PDFs with 
parton flavor a = q,q,g, and Ci/a = X]n=o ^i"a Z'"^)" perturbatively calculable coefficient functions for finding 
a parton i from a parton a, which are given in Ref. fj\. 

The non-perturbative function F^^ in Eq. ([2|) has the following functional form, 



F {b,Mz,XA,XB;b„iax)=exP' 



32 ih" - b: 




(6) 



where the explicit ln(M| dependence was derived by solving CoUins-Soper equation the 52 term is a 

result of adding a general power correction to the renormalization group equation, and the §2 term represents the size 
of intrinsic transverse momentum of active partons f?| . 

For a given choice of bmax and Z° mass, coefficients of two terms proportional to b^ in Eq. ^ can be combined 
together as [27| , 



G, = In 



^"Z"max 



92 + 92: 



(7) 



which sums up both the dynamical and intrinsic power corrections. By requiring the first and second derivatives of W 
to be continuous at 5 = bmax, the parameters: a and gi in Eq. ([5]) can be uniquely fixed leaving only one parameter, G2, 
sensitive to the power correction and other nonperturbative effects. It was found [13] that calculated pr-distribution 
is insensitive to the choice of b^ax if it varies between 0.3 GeV~^ and 0.7 GeV~^. Using g2 = 0.25 ± 0.05 GeV^, and 
52 = 0.01 ± 0.005 GeV^, and b^ax = 0.5 GeY'^, effectively, = 0.324 GeV^ in p+p collisions, theoretical 

predictions using Eq. ([T]) are consistent with all data from Tevatron and the LHC see Fig. [l]in which we see a 
very good agreement between the theoretical calculations and both the CDF ^2§\ and CMS fTo| experimental data 
at = 1.8 TeV and ^/s = 7 TeV, respectively. 
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FIG. 1. Theoretical calculations as in Eq. ([1} for production in p+p (or p+p) collisions are compared with experimental 
data from CDF ;28] and CMS 'If^ at = 1-8 TeV and ^ = 7 TeV, respectively. 



For proton-nucleus collisions, we use EPS09 or nDS nuclear PDFs as parametrized in Refs. [T^ and [T^ to take 
care of the leading twist nuclear effect in PDFs. Although both parametrizations are based on the global fitting of 
experimental data, there is a major difference: EPS09 has an anti-shadowing region for gluon distribution inside a 
large nucleus, while nDS has no such region. We will see immediately the consequence of such a difference in the 
nuclear modification of Z'^ production. At the same time, there could be nuclear size enhanced multiple scattering 
effects. Such nuclear effects can be taken into account through the inclusion of power corrections in the formalism. 



4 



Following the method proposed in Ref. [23|, we take into account these nuclear size enhanced multiple scattering 
effects by letting 32 92^^^^- For Z° production, we have effectively, G^^'' = 0.689 GeV^ in p+Pb collisions at 
LHC. 

In Fig. [21 we plot our predictions for Z'^ production at the LHC energy ^/s — 5 TeV and mid-rapidity y — 
for the planned p+Pb run by evaluating the cross section with resummation in Eq. ([T]) with CTEQ6M NLO parton 
distribution functions and the factorization scale /i — Mt/^ — \/ + We evaluate the F-term at NLO in as 

0. The left panel is for EPS09 nuclear PDFs, while the right panel is for nDS nuclear PDFs. For both panels, the 
upper plots are for the cross section of Z'^ production as a function of its transverse momentum. The black dashed 
curve is the p+p baseline, and the red solid curve is for the minimum bias p+Pb collision. The blue dotted curve is 
exactly the same as the red solid curves except that 52 is not enhanced by A^^^. That is, the nuclear size enhanced 
dynamical power corrections from multiple scattering is absent in blue dotted curve. Lower plots are for the nuclear 
modification factor RpA in min.bias p+Pb collisions. 
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FIG. 2. boson production in both p+p and p+Pb collision at center of mass energy -y/s — 5 TeV and rapidity y = 0. The 
left panel is for EPS09 p^ . while the right panel is for nDS [l^. For both panels, the upper plots are the cross section as a 
function of transverse momentum pr- The black dashed curve is the p+p baseline, and the red solid curve is for the minimum 
bias p+Pb collision. The blue dotted curve is exactly the same as the red solid curves except that 32 is not enhanced by A^^^ . 
That is, the nuclear size enhanced dynamical power corrections from multiple scattering is absent in blue dotted curve. Lower 
panel is for the nuclear modification factor, RpA- 



From Fig. [21 it is clear that the red solid curves are almost indistinguishable from the blue dotted curves. That 
is, the power correction is indeed not important for Z^ boson production at the LHC energy, and therefore, the Z'^ 
production in p+A collisions is an ideal process to probe nuclear modification of parton distribution functions, and 
to test the "isospin" effect at high energies. At the same time, we also find that the cross section of Z^ production at 
^/s = 5 TeV and j/ = is dominated by the gluon initiated subprocesses when pT > 20 GeV. This can be clearly seen 
in Fig. [21 in which we plot the relative contributions (ratios) for Z^ production in p+p collisions from the processes 
which involve gluon distribution (the red solid curve) and those which do not (the blue dashed curve). That is, RpA 
is an excellent observable to measure the nuclear modification of gluon distribution, which is effectively unknown. 

Before we try to analyze and understand the nuclear modification factor RpA of Z^ production in p+Pb collisions. 
Let us first take a look at the difference between EPS09 and nDS nuclear PDFs. In Fig. [31 we plot the nuclear 
modification ratio Rf{x,ii) at factorization scale /i = Mz inside a lead nucleus for both valence u-quark and gluon 
distribution for EPS09 and nDS nuclear PDFs. Even though the valence u-quark distribution of both parametrizations 
is very similar, the gluon distribution is quite different: gluon distribution of EPS09 nuclear PDFs [18\ shadows when 
X < 0.005 while anti-shadows for a sufficiently large range x — (0.005,0.2); gluon distribution of nDS nuclear PDFs 
shadows rather mildly at a; < 0.01 and has almost no nuclear modification for x = (0.01, 0.4). Particularly there is no 
anti-shadowing region for gluon distribution in nDS parametrization. 

Let us now take a look at the nuclear modification factor RpA for production in p+Pb collisions at ^/s = 5 TeV 
and mid-rapidity y = 0. As indicated by the lower panel in Fig. [21 (left), the nuclear modification factor RpA for Z^ 
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FIG. 3. The fractional contributions to production in p+p collisions at -^s — 5 TeV and rapidity y = from partonic 
subprocesses with at least one initial-state gluon (the red solid curve) and those not initiated by gluons (the blue dashed curve). 
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FIG. 4. Ratio of nuclear PDFs over the proton PDFs at scale Q = Mz- Rf ~ fi/A{^, Q"^) /Ji/vjx, Q^)- The red solid curve is 



for gluon, and the blue dashed curve is for the valence M-quark. Left panel is for EPS09 [18| ] and right panel is for nDS [1 



production in EPS09 parametrization is suppressed at low px < 10 GeV, while it is enhanced at high pT, even at 
Pt as large as 70 GeV. The suppression at low px was an immediate result of the shadowed EPS09 nuclear PDFs at 
small X, while the strong enhancement for such a large range of px was a surprise. After a careful examination of 
kinematics, we find that at y = the minimum parton momentum fraction x — Mz/^/s ~ 0.018, which is already in 
the anti-shadowing region of nuclear gluon distribution of EPS09. That is, the clear enhancement of RpA in large px 
region in Fig. [2] (left) is caused by the large anti-shadowing region of nuclear gluon distribution of EPS09. On the 
other hand, for the nDS parametrization, we find that the enhancement of RpA in the large px disappears, as clearly 
seen in the lower panel in Fig. [2] (right). This is due to the absence of (or much smaller) gluon anti-shadowing of 
nuclear PDFs in nDS parametrization. 

In conclusion, the RpA of Z° production is a direct measurement of nuclear gluon distribution since the cross section 
is dominated by gluon initiated subprocess when px > 20 GeV. In particular, the measurement of nuclear modification 
factor RpA of Z'^ production in p-l-Pb collisions at the LHC provides a clean and unique test of the nuclear gluon 
anti-shadowing proposed in EPS09. 
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III. TRANSVERSE MOMENTUM BROADENING OF VECTOR BOSON (J/i/', T, WjZ^) PRODUCTION 

IN PROTON-LEAD COLLISIONS AT THE LHC 



In this section, we will study the transverse momentum broadening of vector boson (J/'0, T, W/Z^) production 
in proton-lead collisions at the LHC. For the hadronic production of inclusive vector bosons, A{pa) + B{pb) — >■ 
F[J/V', T, W/Z^]{q) + X at the LHC, we define the averaged transverse momentum square of the produced vector 
boson as, 



(4X2/) 



AB 



2 2 daAB^V 



dqrp Qrp 



dy dq%, 



dq. 



2 dcTAB^V 

dy dq%. 



We further define transverse momentum broadening in proton-lead collisions as, 

^{qT)pPb{y) = {<iT){y)pPb - (9T)(y)f 



(8) 



(9) 



i.e., the difference in the averaged transverse momentum square of vector boson production in p+Pb and p+p collisions. 
Transverse momentum broadening A{q^)ppi, is directly sensitive to the parton multiple scattering inside a large 
nucleus, as shown in Refs. [23 - l23 |. 

Following the derivation in Refs. [22l. [2^. we have the first nonvanish contribution to the transverse momentum 
broadening of heavy quarkonium production as 



,2 \CEM 



HQ 



m - 1 



^2^1/3 



{Cp + Ca) CFqq + 2 Ca CTgg + Act 



'99 



(10) 



where the superscript "CEM" indicates that the formation of heavy quarkonium from a produced heavy quark pai r 
is evaluated in terms of the Color Evaporation Model (CEM) [2^ , similar result was derived in NRQCD model [l^] . 
The aqq and agg are contributions from quark- antiquark and gluon-gluon subprocess, respectively |22| . 



F, 



QQ-^H 



F, 



QQ^H 



4M,< 



4M'; 



dQ^Yl / Ipq/pix') j dx(j)qlA{x) 

da 



da, 



qq~^QQ 
dQ^ 



dQ / dx'(j)g/p{x') / dx(j>g/A{x) 



dQ^ 



(11) 
(12) 



IS a non- 



where mg is the heavy quark mass, Mq is the mass of the lowest open heavy flavor meson. Fqq^j^ 
perturbative transition probability for a heavy quark pair QQ transforming into a heavy quarkonium H and is 
independent of the color and angular momentum of the heavy quark pair, '^'^''jq^^ and '^'^^jq^P^ represent the total 



cross sections for qq 



and gg — > QQ channels given in |22l |. We reproduce here for convenience. 



dQ^ 

^'^99-^QQ 



iv2 3g2 

1 TTQ!? 



1 



1 + 2^ ) VI 



dQ^ 



1 +7-f 



In 



- 1) ^ 6(iV|-l) V ^ 



l-\/T^ 



(13) 



(14) 



where iVc = 3 is the number of colors, and 7 = AniQ/Q^. 

The AcTgg is a small color suppressed correction to the gluon-gluon subprocess first derived in Ref. [l^l S It is given 

by 



with ^^f££^ 



dQ2 



given by 



dAa 



gg^QQ 



dQ^ 



1 ira' 

nI~T~q^ 



dQ'^<l>g/p{x') / dx<j)g/A{x) 



dAc 



gg- 



dQ^ 



-7-27 



lnii^-(l + 7)v^' 



(15) 



(16) 



^ This piece was missed in the original paper Ref. [22ll . 
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In the region where the gluon-gluon subprocess dominates the heavy quarkoni um p roduction rate, Ugg 
the broadening of heavy quarkonium production could be further simphfied as |22| 



/HQ 



2Ca 



m-1 



> aqq, Aagg, 



(17) 



In Fig. [51 we plot our predictions for the transverse momentum broadening of Drell-Yan type vector boson production 
in p+Pb collisions at the LHC at -y/s = 5 TeV and ?/ = as a function of iVcou, the number of collisions. To get 
the dependence on iVcon, effectively, we replace the A^^^ in Eq. ([TU]) by A^^-^ Ncoiiib) /Ncoiiibmin.his^s) ■ For the p+Pb 
collisions at the LHC, A^coiK^min.bias) ~ 7 from a Glauber model calculation with nucleon-nucleon inelastic cross 
section a™^^ — 70 mb at ^/s — 5 TeV [s^ . For the broadening of heavy quarkonium production, we assume that 
physical heavy quarkonia are formed outside the nucleus, or equivalently, we assume a maximum broadening from 
the final-state multiple scattering. If the color of produced heavy quark pairs is neutralized far within the size of 
the lead nucleus, the rate of multiple scattering of the pairs with the nuclear medium should be reduced, so as the 
broadening, due to the color singlet nature of the pairs. That is, the difference between our predictions and the 
measured broadening of heavy quarkonia at the LHC should provide valuable information on the color nature of 
produced heavy quark pairs and the time scale of color neutralization of the pairs when they transmute into physical 
quarkonia. 

In Fig. [5l we also plot the broadening of W/Z'^ production in p+Pb collisions at the LHC by using the broadening 
formalism derived in Ref. [22]. The main difference in the transverse momentum broadening between the heavy 
quarkonium {J/i/j, T) and W jZ^ production is caused by the following facts: (1) W jZ'^ production has only initial- 
state multiple parton scattering, while heavy quarkonium can receive contributions from both initial-sate and final- 
state multiple parton scattering; and (2) the multiple scattering effect is enhanced for heavy quarkonium production 
because it is dominated by gluon initiated subprocesses. Thus, the dramatic difference in broadening between the 
heavy quarkonium and W jZ^ production in Fig.[5]should be a signature prediction of QCD calculations. The difference 
is sensitive to how a heavy quark pair transmutes itself into a physical quarkonium, and in particular, how the color 
of the pair is neutralized during the hadronization process, which has a direct impact on the final-state multiple 
scattering and its contribution to the amount of the broadening. 
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FIG. 5. Transverse momentum broadening of vector boson production in p+Pb collision is plotted as a function of number of 
collisions Afcoii, we choose LHC kinematics = 5 TeV and mid-rapidity y = 0. Red solid is quarkonium T, red dashed is J/V', 
black solid is W . and black dashed is boson. 



IV. CONCLUSIONS 



In this paper, within the CSS resummation formalism, we first studied the nuclear modification of transverse 
momentum distribution of Z^ production in p+Pb collisions at LHC energy sfs = 5 TeV and mid-rapidity ?/ = 0. We 
found that the production rate of in this kinematic region is dominated by the gluon initiated subprocesses when 
transverse momentum px > 20 GeV. Consequently, the nuclear modification factor RpA is an excellent observable 
to measure the nuclear modification of gluon distribution in a large nucleus, which is practically unknown. Very 
different nuclear gluon distributions were extracted from various QCD global analyses of nuclear PDFs. In this paper, 
we presented predictions for nuclear modification factor RpA of Z'^ production as a function of pT in p+Pb collisions. 
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which will soon be confronted by new experimental data in the upcoming proton-lead run at the LHC. In particular, 
we found that with EPS09, nuclear modification factor is suppressed at low px < 10 GeV, while it is enhanced at 
high pt, even at pT as large as 70 GeV. This clear enhancement of RpA in large pr region is caused by the large 
anti-shadowing region of nuclear giuon distribution of EPS09. The large enhancement at large px actually disappears 
if we use another parametrization of nuclear PDFs (such as uDS) which has much smaller (or no) gluon anti-shadowing 
in the nuclear gluon distribution. Therefore, the measurement of nuclear modification factor RpA of Z° production in 
p-|-Pb collisions at the LHC provides a clean and unique test of the nuclear gluon anti-shadowing proposed in EPS09. 

We also made predictions for the transverse momentum broadening of vector boson {J/ip, T, W/Z'^) production 
in proton- lead collisions at y/s — 5 TeV at the LHC. We found that the transverse momentum broadening for heavy 
quarkonium production is much larger than those for W/Z^ production, due to the difference in color flow of the 
dominated partonic subprocesses for the production of heavy quarkonia and heavy W/Z'^ bosons, and the additional 
final-state parton multiple scattering for heavy quarkonium production. The dramatic difference in broadening be- 
tween the heavy quarkonium and W/Z^ production is a signature prediction of QCD calculations of parton multiple 
scattering in a better-controlled cold nuclear medium. 
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Note added. When this paper was being prepared, a preprint |3l| has just appeared which also studied the nuclear 
modification of Z° boson production at the LHC. Ref. [3l| focused primarily on the small- a; or shadowing region, 
while we study the anti-shadowing region. 
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